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ABSTRACT 

We propose that gamma ray bursts (GRBs) are produced by a shower of heavy blobs running into 
circumstellar material at highly relativistic speeds. The gamma ray emission is produced in the shocks 
these bullets drive into the surrounding medium. The short term variability seen in GRBs is set by 
the slowing-down time of the bullets while the overall duration of the burst is set by the lifetime of the 
central engine. A requirement of this model is that the ambient medium be dense, consistent with a 
strong stellar wind. The efficiency of the burst can be relatively high. 
Subject headings: gamma rays: bursts 



1. INTRODUCTION 

Gamma ray bursts (GRBs) show variability over a large 
range in time scales — from millisecond spikes in BATSE 
light curves (Bhat et al. 1992) to months in the associated 
afterglows (Bloom et al. 1999, Frail et al. 1997). The dis- 
covery of afterglow redshifts that place GRBs at cosmolog- 
ical distances implies isotropic energies of E mo <; 10 53 ergs 
for the average burst in gamma rays alone, released over 
an observed time of order T ^ s ~ 10 s. Since GRB spectra 
are believed to be optically thin, the only viable explana- 
tion is a highly relativistic outflow that releases its energy 
at sufficient distances from the central engine for the opti- 
cal depth and the compactness of the plasma to be small, 
with required bulk Lorentz factors of T ^ 100. 

In early models (Meszaros & Rees 1992, Katz 1994a, 
Sari & Piran 1995) the gamma rays were produced by the 
shock such an outflow drives into the interstellar medium 
(thus called 'external shock models'). While providing a 
way to produce the required energy in gamma rays over 
the observed GRB durations, these models have not been 
very successful at explaining the observed short term vari- 
ability: in order to see strong variability, the surrounding 
medium must be very clumpy (e.g., Shaviv & Dar 1995, 
Fenimore et al. 1996). It has been argued that this pro- 
cess might lead to very low efficiencies at converting kinetic 
into thermal energy (Sari & Piran 1997). However, a quan- 
titative investigation casts some doubt on this argument 
(Dermer & Mitman 1999). 

These difficulties of the external shock model prompted 
Sari & Piran (1997) to postulate that the gamma ray 
emission must instead be produced by the internal shock 
scenario (Narayan, Paczyhski, & Piran 1992, Rees & 
Meszaros 1994). In this picture, the duration of the burst 
T is set by the time scale over which the central engine op- 
erates, while the substructure in the bursts on time scales 
r is produced by inhomogeneities in the outflow. These 
inhomogeneities are assumed to travel at different bulk 



Lorentz factors T. Upon running into each other, these 
shells of material get shocked and release some of their ki- 
netic energy in the form of gamma rays. This picture has 
become the paradigm in GRB physics. However, recent es- 
timates of the energy conversion efficiency rj indicate that 
at best a few percent of the bulk kinetic energy carried by 
the outflow can be converted into gamma rays in internal 
shocks, which leads to uncomfortably high requirements on 
GRB energies (Panaitescu, Spada, & Meszaros 1999, Ku- 
mar 1999; larger efficiencies can be achieved if you assume 
large dispersion in T, Katz 1997). A non-spherical geome- 
try can reduce the required energy, however, a very small 
opening angle of the outflow implies a high rate of un- 
observed GRBs, which is hard to reconcile with the num- 
ber of observed supernovae (which are believed to produce 
GRB precursors - either compact objects or hypernovae). 
See Piran (1999) for a detailed review of the subject. 

Fenimore et al. (1996, 1999, see also Woods & Loeb 
1995) recently suggested that an external shock scenario 
could give rise to the observed variability if the spherical 
symmetry of the outflow were broken, still in the context 
of what Sari & Piran (1997) call a 'Type I' model, i.e., 
the burst duration is set by the slowing-down time of the 
ejecta. However, the observed temporal constancy of the 
pulse width in individual spikes of GRB990123 seems to 
rule out such a model (Fenimore, Ramirez-Ruiz, & Wu 
1999. This would also argue against the scenario suggested 
by Dermer & Mitman, 1999.) We propose a different way 
by which substructure in the outflow can produce a GRB, 
also via the interaction with the external medium, but in 
a 'Type IF scenario, i.e., the duration of the burst is set 
by the lifetime of the central engine (see also Chiang & 
Dermer 1999, Blackman, Yi, & Field 1996). 

In this model, the outflow itself is very clumpy, with 
most of the energy concentrated in small blobs, which are 
sprayed out with high T over a small opening angle. These 
bullets then slam into the surrounding medium (not unlike 
a meteor shower or a shotgun blast), where they release 
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their kinetic energy and produce gamma ray emission via 
external shocks, as described in §2. Note that, in this pa- 
per, we assume the existence of bullets; their possible ori- 
gin is a topic for further research (see §3 for comments). 
Section 3 contains a discussion on radiative efficiencies, 
simulated burst light curves, and afterglows, and §4 con- 
tains a brief summary. 

2. THE GAMMA RAY BURST MODEL 

BATSE light curves of most GRBs exhibit very spiky 
substructure (Fig. 1). In the context of our model, these 
spikes are produced by individual bullets of cold ejecta 
slamming into the surrounding medium. As we will show, 
a distribution of masses and/or Lorentz factors of these 
bullets can reproduce the observed signatures of GRBs 
reasonably well. 
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Fig. 1— BATSE light curves for GRB920627 (panel a) and 
GRB980329 (panel b) and two synthetic light curves. These curves 
were calculated for a burst duration of 15 s, a mass distribution of 
N(Mb) oc Mb -1 / 3 , an average slowing down time of r = 0.01 s, and 
^100 = 1 (panel c), A?ioo = 100 (panel d). The insert in panel c 
shows a template light curve for a single shot. 

In the following we assume that the central engine of 
the burst releases a number of bullets N' distributed over 
a fan of opening angle of 6 <~ 10° with Lorentz factor T = 
1000 and released over a time period of T ~ 10 — 100 s. 
Each bullet is assumed to freely expand sideways, with a 
sideways velocity of v± = a/T = 10~ 2 a_2/r <C 1/T, mea- 
sured in the observer's frame (a is the sideways velocity in 
the comoving frame). The assumption that a< 1 implies 
that the internal expansion speed is very sub-relativistic. 

Since 6 » 1/T, we only see a fraction of the total re- 
leased energy, -E b s = 10 47 • E53/T3 2 ergs (where -E53 is the 
inferred isotropic energy in units of 10 53 ergs) and an ob- 
served number of bullets N = 100 • N loa = A'/(6l 2 r 2 ). It 
is essential in our model that the covering fraction be less 
than unity (otherwise it would turn into an internal shock 
model), thus we require a < 1/y/N. 

A bullet of mass M b = 5.5 x lO- 13 ^ 5 3M /7V loo r 3 3 will 
have converted half of its kinetic energy into internal en- 
ergy (which can subsequently be radiated away as gamma 
rays, see below) after it has swept up or ploughed through 
a column of interstellar gas of mass M s ~ Mb/T. If the 
duration of the observed spike is t s = 10 _3 t_3 s, then the 
length over which the material is swept up must be of the 



order of AR 

is n am b = 10 8 n§ cm -3 , the required Lorentz factor is 
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This model fails for low ambient densities, as has al- 
ready been discussed in the literature (e.g., Sari & Piran 
1997). However, if the surrounding medium is very dense, 
n ~ 10 8 cm~ 3 , Lorentz factors of F <~ 1000 can explain the 
observed short term variability. The immediate conclu- 
sion is that in this model GRBs are not caused by merg- 
ers of naked compact objects. Rather, the required high 
ambient densities tie this model to the hypernova picture 
(Paczyhski 1998, Woosley 1993), which predicts that the 
material surrounding the blast is dense because of the prc- 
hypernova stellar wind. (Dense circum-GRB matter was 
also suggested in a different context by Katz, 1994b.) 

The fate of the outer layers of a hypernova precursor 
is unknown. If the bullets have to travel through a sig- 
nificant fraction of the star's mantle (which is optically 
thick and thus not the site where the gamma rays are pro- 
duced), their opening angle must be extremely narrow: 
a < 5 x 10~ 4 S 5 3/(v^Vi~mP3-Mm), where M M is the mass of 
the mantle in units of Mq. Since hypernovae are believed 
to originate from rapidly rotating massive stars collaps- 
ing into a compact object, it is possible that the material 
along the rotation axis has collapsed before the GRB, in 
which case the bullets would travel freely until they hit 
the circumstcllar material. 

Similarly, little is known about the conditions of the pre- 
hypernova circumstellar material other than that it must 
be dense. Massive stars are known to have strong winds 
with mass loss rates of M <~ 10 -6 — 10 _4 M Q yr _1 and 
wind velocities from v w ~ 20kms _1 (red supergiant) to 
w w ~ 1000 km s -1 (blue supergiant). These winds must 
still be present after the star collapses. In the following, 
we will assume that the GRB is produced in this leftover 
wind. The density profile in the ambient matter, then, 
roughly goes as n cx r~ 2 outside some radius Rq. It is nat- 
ural to assume that Rq is of the order of the stellar radius, 
i?o ~ 10 12 cm for a blue supergiant and Rq <~ 10 14 cm for 
a red supergiant. As a conservative estimate we assume 
that the sphere inside Rq is evacuated. 

If we parameterize the density as 



n amh (r > Rq) 



1.5 x 10 36 cm- 3 M_ 6 
r 2 v 20 
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where M-e is the mass loss rate in units of 10 _6 M Q yr _1 , r 
is in cm, and W20 is the wind velocity in units of 20kms~ 1 , 
the observed slowing-down time scale is given by 



r = 0.05 s- 
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independent of i? - Thus t can be of the order of a 
few milliseconds for both red supergiant and blue super- 
giant winds if T ~ 1000. However, the observed time 
scale could conceivably be longer than this. The an- 
gular smearing time scale r ang is defined as the spread 
in light travel time to the observer across the emitting 
surface. For a bullet at a viewing angle of 1/r this 
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gives T ang ~ a(Ro + AR)/T 2 , which is longer than r if 
Ro^ AR/a = 1.5 x 10 17 cm - E 53 V2o/{r 3 2 a-2 3 N wo M r6 ). 
This is only of concern for very dense red supergiant winds, 
and only if the region interior to Rq is evacuated. 

In Fig. 2 we show various limits on n am b and T for a 
fixed opening angle of a -2 = 1' 

a) Each bullet is expected to plough through undis- 
turbed medium. Thus, the covering fraction Na 2 
of all the bullets together must be smaller than 1. 
For a slowing down time of r_3 = 1, this gives the 
dashed line in the plot, to the left of which the cov- 
ering fraction is larger than unity. 

b) The material between the location where the bul- 
lets release their energy and the observer must be 
optically thin. For r_3 = 1 and for Rq = (the 
most conservative limits) this constraint produces 
the dash-dotted line to the left of which the opti- 
cal depth is larger than unity. 

c) The forward shock must be radiative (see §||). This 
constraint is shown as a light grey region inside of 
which the shock is not radiative. 

d) The angular smearing time r ang must be smaller than 
the observed stopping time r. This limit is shown as 
a dotted line for Rq = 10 14 cm. To the right of this 
line, the smearing time is longer than the observed 
slowing down time. 

The hatched region in the plot shows how the allowed re- 
gion of parameter space opens up if we relax the time scale 
requirements to r_3 = 10. 




r 

Fig. 2. — Constraints on the ambient density and T for £53 = 1, 
Niqq = 1, and a = 0.01. a) covering factor larger than unity: left 
of dashed line for r_3 = 1, dark grey area for t_3 = 10. b) optical 
depth between bullet and observer larger than unity: left of dash- 
dotted line for t_3 = 1, medium grey area for r_3 = 10. c) forward 
shock not radiative: light grey area (we assumed eb = !)■ d) an- 
gular smearing time longer than slowing down time: right of dotted 
line. The hatched area shows how the allowed region in parameter 
space expands if we relax the requirement on r from t_3 = 1 to 
T-3 = 10. 



It is worth noting that this model makes an exception 
to the rule that external shocks cannot produce 'Type II' 
behavior (Sari & Piran 1997). This is for two reasons: 
First, the opening angle of the ejecta is so small that the 
angular smearing time is short compared to r. Second, 
the ambient density is so high that the observed slowing- 
down time is t ~ lCP 3 s. As a result, the total duration 
of the burst is determined by the time the central engine 
operates, while the short term variability is determined by 
the mass of the bullets and the statistical properties of the 
outflow. This is an important difference from the internal 
shock model, where the variability timescale is set by the 
intrinsic time scale of the central engine (e.g., the orbital 
time in a merger scenario). 

3. DISCUSSION 

In order for the efficiency of the burst to be reasonable, 
most of the internal energy produced in the shock must 
be radiated away immediately (this requirement holds 
for all GRB models). Electron synchrotron radiation is 
the only mechanism remotely efficient enough to produce 
the gamma rays. While the efficiency also depends on 
the transfer of energy from protons to electrons, we as- 
sume here that this process is efficient. Since we know 
the observed peak frequency of the gamma rays (of or- 
der 500 keV, Piran 1999), we can then estimate the radia- 
tive efficiency under the assumption that the gamma rays 
are produced by electron synchrotron radiation. It is usu- 
ally assumed that the magnetic field in the shocked gas 
is in equipartition with the energy density in relativistic 
particles. We therefore parameterize the magnetic field 
strength as Ub = eB^B,cq, where i7seq is the equipartition 
magnetic field energy density. 

Since the shocked wind material is likely flowing around 
the bullets at close to the speed of light (like a cocoon sur- 
rounding a radio jet), for efficient cooling we require that 
the cooling time scale in the comoving frame be smaller 
than the light travel time across the surface of the bullet 
a(i?o + AR)/Tc. If that were not the case, the material 
pushed aside by the bullet would cool adiabatically rather 
than radiatively. Independent of Rq, this translates to the 

condition T < 1430 • [M_ 6 a_ 2 2 £B Mo] 1/4 \/^/ 500 kcV ( to 
the left of the light grey area in Fig. 2), where E v is the 
observed peak energy of the gamma rays. This is not a 
strict condition, however, since we do not know what the 
efficiency of the GRB is. 

Since our model assumes a central engine at work (es- 
sentially a black box shooting out bullets at a rate R(i)), 
any distribution of spikes could be reproduced, since R(t) 
is arbitrary. It is, however, surprisingly simple to repro- 
duce the main features seen in different burst profiles by 
varying only a few parameters in our model. For simplic- 
ity, we assume that all the bullets have the same initial 
r and the same sideways expansion rate, i.e., constant a. 
We are left with two parameters — the number of bullets 
-/Vioo and the average slowing-down time (eq. ||) — and 
two unknown functions: the mass distribution of the bul- 
lets iV(Mb) and the rate at which they are released R(t). 
We assume that N(Mb) oc Afb 1 ^ 3 , chosen to give the 
observed power spectrum of P{t) oc t -1 / 3 (Beloborodov, 
Stern, & Svensson 1999). For R{t) we assume (for lack of 
better knowledge) that the bullets are released randomly 
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over a time interval of 15 sec. 

To produce synthetic GRB lightcurves, we calculated 
the time dependence of T and the associated dissipation 
rate. Assuming the bullets are radiating efficiently and 
correcting for Dopplcr boosting and frequency shifts, we 
then computed the composite lightcurve for each bullet. 
We have plotted two simulated light curves in Fig. 1 for 
t_ 3 = 10 and N = 100 (panel c) and N = 10 4 (panel 
d). It seems that simply by varying Aioo and t_3 we can 
produce a wide range in light curve shapes. More complex 
features (like the gaps seen in panel a) must be related 
to the activity of the central engine and cannot be repro- 
duced by a random spike rate as assumed above. We have 
also plotted the light curve produced by the deceleration 
of a single bullet along the line of sight (Fig. lc, insert). 
Note that this profile is very similar to a true FRED (fast 
rise, exponential decay) profile. The rise is instantaneous 
and the decay follows a steep power-law (to first order). 
While indicative, these calculations are still rather crude 
and simplistic. A more careful analysis of shotgun GRB 
lightcurves should be carried out in the near future. 

Afterglows are an important test for any GRB model. 
How can we understand an afterglow in the context of the 
shotgun model? It is not immediately obvious why our 
model should produce an afterglow at all. This is because 
the bullets are assumed to spread sideways. As mentioned 
above, the bullets will have lost half of their kinetic en- 
ergy at AR, where they have swept up 1 /T of their own 
mass. If we simply followed the dynamics of an individual 
bullet further in time, it would lose the rest of its energy 
exponentially quickly (Rhoads 1997). This is because the 
sideways velocity in the lab frame goes as 1/r, so that 
when r starts decreasing, the sideways velocity increases, 
which in turn increases the cross-sectional area of the bul- 
let. As a result, the bullet can sweep up more mass, which 
leads to a run-away process. This would imply that the 
ejecta would come to a complete stop not far away from 
A_R and the observed afterglow would last less than a day. 

However, there are many bullets traveling together. As 
they expand, they start increasing the covering factor of 
the blast. Once the collection of bullets reaches unit cov- 
ering factor, they stop slowing down exponentially since 
further sideways expansion does not lead to an increase 
in swept up mass. As a result, the bullets start travel- 
ing collectively, resembling a collimated blast wave with 



opening angle 9 rather than a meteor shower. The only 
possible difference between our model and the standard af- 
terglow models is that in our model, the external density 
follows a power-law behavior instead of being constant, 
which has been discussed by Dai & Lu (1998). Since the 
opening angle of the merged blast is much larger than 
the opening angle of the individual bullets, the sideways 
expansion does not affect the dynamics until much later, 
when the blast has spread by an angle of order 9 (Rhoads 
1997). This transition from constant opening angle to 
rapid sideways expansion has been used in other models to 
explain the temporal break seen in the afterglow lightcurve 
of GRB990510 (Stanek, et al. 1999, Harrison et al. 1999, 
Sari, Piran, & Halpern 1999) 

It should be noted that we have not attempted to 
present a specific model for the production of bullets in 
this paper. We are confident that these bullets can be 
produced in GRB outflows and will concentrate on their 
production in future work. 

4. CONCLUSIONS 

We have shown that GRBs can be produced by a shower 
of cold, heavy bullets shot at bulk Lorentz factor Y <~ 1000 
into a dense medium. The required densities are consistent 
with a stellar wind from either a blue or red supergiant. 
This tics our model directly to the hypernova scenario. 
The gamma rays are produced by the shocks these bul- 
lets drive into the ambient gas. The total duration of the 
burst is then determined by the time the central engine 
operates rather than the slowing-down time of the bullets, 
while the latter produces the short-term variability seen in 
many bursts. After the gamma ray phase (when the bul- 
lets have lost half of their kinetic energy to radiation) the 
blast waves of the individual bullets merge into a single 
collimated shock front, which produces a standard after- 
glow in a declining external density profile. We have made 
no attempt to explain how a central engine might produce 
such a shower of bullets. This will be the subject of future 
work. 
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